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Abstract. We develop the theory of a metamaterial composed of an array of 
discrete quantum absorbers inside a one-dimensional waveguide that implements 
a high-efficiency microwave photon detector. A basic design consists of a few 
metastable superconducting nanocircuits spread inside and coupled to a one- 
dimensional waveguide in a circuit QED setup. The arrival of a propagating quantum 
microwave field induces an irreversible change in the population of the internal levels 
of the absorbers, due to a selective absorption of photon excitations. This design is 
studied using a formal but simple quantum field theory, which allows us to evaluate the 
single-photon absorption efficiency for one and many absorber setups. As an example, 
we consider a particular design that combines a coplanar coaxial waveguide with 
superconducting phase qubits, a natural but not exclusive playground for experimental 
implementations. This work and a possible experimental realization may stimulate 
the possible arrival of " all-optical" quantum information processing with propagating 
quantum microwaves, where a microwave photodetector could play a key role. 
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1. Introduction 

In a recent work pQ we suggested a possible implementation of a photon detector that 
may also work as a photon counter in the microwave regime. Our proposal builds on 
previous advances in the field of quantum circuits in two fronts. One is the development 
of artificial atoms and qubits [21 El HI [5] for quantum computation and quantum 
information processing, using quantized charge [6j [7J [U [9], flux [TQl [HI [T21 [13], or 
phase [H], [15], [TU [T7] degrees of freedom. The other front is the efficient coupling of 
these elements to microwave guides and cavities conforming the emergent field of circuit 
quantum electrodynamics (QED) [HI [191 EO] . Without neglecting important advances 
in intracavity field physics in circuit QED, as we will continue to illustrate below, here 
we are interested in the physics of propagating quantum microwaves. 

As we argued before pQ , in order to fully unleash the power of quantum correlations 
in propagating microwave photonic fields, as may be generated by circuit QED 
setups, the implementation of efficient photon detectors and counters would be mostly 
welcomed. The existence of these detectors is implied by almost any sophisticated 
quantum protocol involving optical photons, be in their coherent interaction with 
matter [21J or purely all-optical devices [22]- It ranges from the characterization 
and reconstruction of nonclassical states of propagating light by quantum homodyne 
tomography [23] to high-fidelity electron-shelving atomic qubit readout [24|. Both 
examples coming from quantum optics have shown to be influential in the novel field 
of circuit QED, with the first theoretical [25J and experimental efforts [26J to measure 
relevant observables of propagating microwaves, and a recent proposal of mesoscopic 
shelving qubit readout [27]. In spite of these efforts, it will be very hard to overcome 
the necessity of photon detectors and counters when the emerging field of quantum 
microwaves will want to deal with local and remote interqubit/intercavity quantum 
communication, implementations of quantum cryptography, and other key advanced 
quantum information protocols [2H [22] . 

It should be thus no wonder that photon detection and counting become soon 
a central topic in the field of quantum circuits, where superconducting circuits 
interact with intracavity and propagating quantum microwaves. So far we have seen 
the exchange of individual photons between superconducting qubits and quantum 
resonators [2H [29], [30], the resolution of photon number states in a superconducting 
circuit |31j, the generation of propagating single photons [32], the first theoretical 
efforts for detecting travelling photons J33[ [M], and the nonlinear effects that arise 
from the presence of a qubit in a resonator [351 136]. We envision a rich dialogue between 
intracavity and intercavity physics in the microwave domain, see for example [37J,|38] 39J, 
where matter and photonic qubits exchange quantum information in properly activated 
quantum networks for the sake of quantum information processing. 

All efforts towards the implementation of a photodetector for propagating 
microwaves in circuit QED face a number of challenges, many of which are related 
to the specific nature of quantum circuits pp. These are: i) Available cryogenic 
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Figure 1. Our photodetector scheme consists of three- level absorbers (a) embedded 
in a coplanar coaxial microwave guide (long gray stripes). The absorbers (squares) 
can either be grouped in clusters that are much smaller than a wavelength (b) or 
regularly spaced (c). As the waveguide provides an effectively one-dimensional setup, 
the transverse position of the absorbers affects mildly the coupling strength. 

linear amplifiers are unable to resolve the few photon regime, ii) Free-space cross- 
section between microwave fields and matter qubits are known to be small, iii) The 
use of cavities to enhance the coupling introduces additional problems, such as the 
frequency mode matching and the compromise between high-Q and high reflectivity, 
iv) The impossibility of performing continuous measurement without backaction |33j, 
which leads to the problem of synchronizing the detection process with the arrival 
of the measured field. In a wide sense, the photodetection device has to be passive, 
being activated irreversibly by the arriving microwave signal. Otherwise, the advanced 
information that a photon is approaching turns itself into a photodetection device. 

Our proposal for a photon detector consists of a very simple setup, a microwave 
guide plus a number of superconducting circuits that absorb photons [Fig. [Q, and is 
able to circumvent at once the problems describe above. Instead of unitary evolutions, 
we make use of an irreversible process which maps an excitation of the travelling 
electromagnetic field (a photon) into an excitation of a localized quantum circuit. By 
separating this encoding process from the later readout of this information, we avoid 
the backaction problem coming from continuous measurement. Furthermore, we can 
compensate different limitations — weak coupling, low efficiency of absorption, photon 
bandwidth — using no more than a few absorbing elements [Fig. [lb], where collective 
effects enhance the detection efficiency. 

In this work we develop in great detail the theory underlying our proposal for 
high-efficiency phodetection pQ. In Sec. [2} we develop an abstract model that consists 
on a one-dimensional wave guide that transports photons and a number of three- 
level quantum systems that may absorb those photons. We will solve analytically the 
evolution of an incoming wavepacket, studying the time evolution of the full system with 
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one or more detecting elements. For a single photon we can compute the absorption 
efficiency and how it depends on various parameters: detuning of photons from the 
absorber transition, irreversible decay rate of absorbers, spatial distribution of detectors 
and temporal profile of photon. In Sec. [31 we go further with this theoretical model by 
applying it to a particular design where the waveguide is a coplanar coaxial microwave 
guide, built with superconducting materials, and the absorbing elements are built from 
superconducting phase qubits. The parameters of the abstract model are matched to the 
experimental frequencies, impedances and capacitances of the chosen setup, verifying 
the experimental viability of our model. For completeness, in Sec. HI we write down 
the quantum field theory for the microwave guide as it is used and needed in Sec. [2] 
and Sec. [31 Finally, in Sec. [5j we summarize our results and discuss the experimental 
challenges and future scopes. 



2. Abstract model and design 

In this section we present in an abstract manner our photodetector design: a one- 
dimensional waveguide transporting photons that can be absorbed by one or more 
bistable elements that, for the sake of simplicity, will be called indistinctively absorbers 
or qubits. We will introduce a general mathematical model that uses two fields to 
describe the quantized waves (photons) and three-level systems for the absorbers. We 
will study the consequences of this model, including how the joint device may act as 
a photon detector or photon counting device. This mathematical formalism will serve 
as a foundation for the developments in Sec. El where we will specify the microscopic 
theory behind the waveguide and the photon absorbers. 



2.1. Master equation 

Following the standard quantum optical master equation formalism, we describe the 
joint state of the propagating photons and the absorbing elements using a density 
matrix, p. This is a Hermitian operator that contains information about the probability 
distributions of all observables in the system. To the lowest order of approximation the 
density matrix will evolve according to the master equation 

m 

Here, the conservative part of the evolution is ruled by the Hamiltonian 

H = //photon + -f/qubit + Hint, (2) 

which contains terms describing the propagation of microwave photons in the guide, 
the absorbing elements and the interaction between matter (absorbers) and radiation 
(field), respectively. The first operator contains two propagating fields, ip r and ipi moving 
rightwards and leftwards with group velocity v 



photon 



ipl(—ihvd x )'ip r + ipj(+ihvd x )ipi dx. (3) 
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The second part models one or more discrete quantum elements that we place close to 
the transmission line. These elements are analogous to the qubits in quantum computing 
and circuit-QED setups, and will play the role of absorbers, or qubits, enjoying at least 
three energy levels. The first two, |0) and |1), are metastable and separated by an 
energy hu close to the frequency of the incoming photons 

TV 

H quhit = J2^\i)AM- (4) 

i=i 

Then, there is the interaction between the electromagnetic field and our qubits. We 
model it with a delta-potential which induces transitions between the qubit states at 
the same time it steals or deposits photons in the wave guide 

N 

H mt = E / VS ( X - + i)i{x)\ (0| dx + H.c. (5) 

1=1 J 

Finally we have included a Liouvillian operator C which models the decay of the 
absorbing elements from the metastable state |1) to a third state, \g) , and which 
constitutes the detection process itself. A general second order Markovian model for 
the decay operator reads 

N T-l 

£P = E 2 P \9)t I 1 )* (9\ ~ \lh (1\P-P <1|] • (6) 
i=i 

Note that if we start with a decoupled qubit (V = 0) in state |1), the population of this 
state is depleted at a rate T 

p(t) = e- rt \l)(l\ + .... (7) 
2.2. Non-Hermitian solution 

Let us consider the simple case of one qubit or absorber. The master equation (pQ) can 
be written in a more convenient form 

j t P = Ap + pA^ + T\g)(l\p\l)(g\, (8) 

where we have introduced a non-Hermitian operator 

A = -ltf-E|l)<l|^t. (9) 

The master equation can now be manipulated formally using the "interaction" picture 
p(t) = e At a(t)e AH } (10) 
with the following equation for cr(t), 

±a = re- M \g)(l\e At o-e AU \l)(g\e- AU . (11) 



Using the relation e At \g) = e At * \g) = \g) we obtain 

^-a = T\g)(l\e M ae AH \l)(g\. (12) 



dt 
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This equation can be solved by splitting into a part that already decayed, a gg oc \g) (g\ 
and all other contributions o~j_ = 1 — o gg 



^ = r|^)(^|x5:(l|e A V ± e Att |l), (13) 

a, (3 



l« 1= 0. (14) 
Using the initial condition 

<r a {0) = 0, o-x(O) = p(0) = |*(0)> (¥(0)| , (15) 
we obtain the formal solution 

p(t) = Mt)) + l*> fol /*r(l|*(r)) (¥(r)|l) dr, (16) 

where we have introduced state | ¥(£)), evolving from initial state |¥(0)) according to 
the lossy Schrodinger equation 

j t \*(t)) = AMt)) = ~HMt)). (17) 
Note that the probability of the qubit irreversibly capturing a photon is given by 

P(p-+l-+g)=Tr(p\g)(g\) (18) 
= 1-Tr[|*(0> (*(f)|] 

= 1-||¥|| 2 . (19) 
2.3. Model equations for a single absorber 

Following the previous reasoning, our goal is to solve the lossy Schrodinger equation 

i±m = H\9), (20) 
ruled by the non-Hermitian operator 

H = Hphotcm + H qubit + H int — X! "o - l"0t (^l • (^"0 

i ^ 

In order to describe the detection process, our initial condition |¥(0)) will contain the 
incoming photons plus all absorbers in the ground state |0). The probability that an 
absorber captures a photon and changes from state |0) to \g) is given by the loss of 
norm shown in Eq. (fl9~]l . 

Let us consider an incident photon coming from the left with energy E = h\k\v. 
The state of the system will be given by [ID] 

I*) = / [Zr(x,t)^t(x)+ti(x,t)ip}(x)] \vac,0) + 

+ ei \vac, 1) . (22) 

Here ipli(x) \vac, 0) is the state of a photon created at position x and moving either to 
the right or to the left, while the absorber is in metastable state |0). Also, £ r (x,t) and 
£i(x, t) represent the wavefunction of a single photon moving to the right and to the 
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Figure 2. An incident photon, moving rightwards, interacts with an absorbing 
element. Out of the original amplitude of the field, a, a component is transmitted, a', 
another component is reflected, b, and finally with some probability, |c| 2 , the system 
absorbs a photon and changes state. 



left, respectively, while \vac, 1) is a state with no photons and the absorber excited to 
unstable level |1) . Note that thanks to the relation ffl6|) we do not need to explicitely 
include the population of state \g) . We only have to solve the Schrodinger equation 
(I20p using a boundary condition that represents a photon coming from the left and an 
inactive absorber 

£ r (x > 0, t ) = ^(x, t ) = 0, ei(t ) = 0, (23) 

and compute the evolution of the photon amplitude, £ r> i(x, t), the excited state 
population e\(t) and the resulting photon absorption probability ( TT9j) . 

After decomposing the wave equation into the left (x < 0) and right (x > 0) halves 
of space, and replacing the potential 8{x — x±) with an appropriate boundary condition 
at X\ = 0, one obtains the Schrodinger equation for the absorber 

1 k,(0 + )+^(0-)+6(0 + )+£;((r 



(24) 



(25) 



id t ei 

All L 

+ (oj — zr)ei 
and four equations for the photon, 

id t £ r (x,t) = - ivd x £ r (x,t), x ^ 0, 
idS{x,t) = +ivd x £i(x,t), x ^ 0, 
= -ihv[£ r {0 + ,t) -£ r (0- *)] +Ve 1 {t), 

= +ihvUo + ,t)-Uo~,t)} + v ei (t). 

We introduce new variables a, b, a' and b' describing the amplitude of the fields on both 
sides of an absorber [Fig. [2] , 

a(t) =£r((T,t), a'(t) = e(0 + ,t), 
b(t)=£ l (p-,t), b'(t)=U0 + ,t). 

Two of these variables can be solved from the initial conditions 

a{t) =£ r (0 + ,t)=U-v(t-t ),t ) (27) 
b\t) = 0, 



(26) 



Photodetecting propagating microwaves in circuit QED 



8 



id t ei = ( lu — iT — i^~ 1 e i + ~ a (^)' (^9) 



while the rest are extracted from the boundary conditions 

a\t) = a{t) + ^ei(t), &(*) = ^ei(t). (28) 

We have expressed all unknowns in terms of e\(t). The population of the excited state 
now satisfies 

V 2 \ V 

Wv-) ei + T 

which contains the original frequency, u, and an imaginary component for the decay 
of the three-level system, T, enhanced by the interaction with the transmission line 
V 2 /h 2 v. Solving for e±, we find 

ei (t)= -{Y- f e- f{t - T) a(r)dT, (30) 
n Jo 

V 2 

r = r + -g- + icu, 

n v 



in terms of the incoming beam (J2] 

Using the expression of exit) one may compute the reflected and transmitted 
components of a microwave beam, to find that the total intensity is not conserved. 
Rather, we have an absorbed component that is "stolen" by the three-level system to 
perform a transition from the |0) to \g) states via |1) . The balance equations are thus 
[Fig.! 

|a| 2 + |6'| 2 = |a'| 2 + |6| 2 + |c| 2 . (31) 

If we consider a photon with a finite duration, T, the detection efficiency will be defined 
as the fraction of the wavepacket that was absorbed 

Io\c(t)\ 2 dt 

ff[\a{t)\* + \V(t)\ 2 ]dt 
For the purposes of photodetection and photon counting we would like this value to 
reach the maximum a = 1. 



2.4- Stationary solutions for one qubit 

To understand the photon absorption process we will study a continuous monochromatic 
beam which is slowly switched on, 

a(t) = exp[(-iw Q + <r)t] t<0,a>0. (33) 

Taking the limits to — > — oo and a —>■ + in this precise order, we obtain 

y /•* 



V ft 

e x — — i— I exp[— r(t - r) + (— ioo + a)r 
h J-oo 

v/h 



e~ iuJot . (34) 

r + p- + i(u - ujq) 



h v 

We will introduce a unique parameter to describe the photodetection process 

hv 

V 



7=— n[T + i(u-Uo)), (35) 
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which includes both a renormalization of the decay rate and a small imaginary 
component associated to the detuning. With this parameter the solution reads 

e 'W = - 'v TT+™' (36) 



«'(() 



i- 1 



1 + 7. 
1 



e~ iuJot , (37) 



bit) = — e- itJot . (38) 

1 + 7 

If we work in the perfectly tuned regime, u = uo, the decay rate 7 becomes real and 
the absorption rate is 

a = - 2Re(a*b) - 2\b\ 2 (39) 
2 2 

= I + 7" ~ (I + 7) 2 (40) 

= ^ = 26(1-6). (41) 

This value achieves a maximum of 50% efficiency or a = 1/2 at the values b = 1/2, 
7 = 1. We think that the limit of 50% in the photodetection efficiency is fundamental 
and related to the Zeno effect, expressing the balance of quantum information between, 
see Fig. la, the reversible absorption of the photon in the first (left) transition channel 
and the irreversible absorption in the second (rigth) one. 

2.5. Transfer matrix 

We can derive the long wavepacket or quasi-stationary solution in a slightly different 
manner. Note that for infinitely long wavepackets the population of the excited state is 
determined by the fields on both sides 

ei = t Vtt ^ [a + a' + b + b'] (42) 

fr V r / 7 7/1 

[a + a' + b + b'} . 



iV'j 

With this the boundary conditions in Eq. ( |28|) transform into a set of equations that 
only involves the incoming and outgoing fields, 

= a'-a+ — (a + a +6 + 6'), (43) 
27 

= 6'-6-— (a + a' + 6 + 6'). (44) 
27 

In terms of the matrix and vectors 



we can write 

(i+ .... 

27 J \ 27 



1 + — A j a/ = ( 1 - — A ] x. (46) 
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Multiplying everything by (1 — A/2~/) and using A 2 = gives x' = (1 — A/j) x. This 
amounts to a relation between the fields to the right and to the left of the absorber 



a' \ „ / a 



\vj T UJ' (47) 

given by the scattering matrix 

r=( x ^ ^J-i-U (48) 

Note that the value 7 introduced before (155]) only depends on the properties of the qubit, 
the interaction between the qubit and the line and the group velocity. We have thus 
elaborated a rather compact and easily generalizable scheme for studying the scattering 
and absorption of photons by the individual absorbers. 

2. 6. General absorption formula 

Using the transfer matrix formalism we can compute the fraction of absorbed photons. 
The only requirement to obtain a general formula is that the scattering matrix remains 
the same under inversion 

CMS)- 

This is reasonable: from Fig. [2] we can conclude that b and a' have the same role, just 
like b' and a. Using the relations b = T ia' and a = T n a', the absorbed fraction becomes 

« = 1 - 1 J = 1 - pp-faCl + l^oil 2 )- (50) 
I "I I ii I 

This is consistent with the single absorber case 

T n = T 01 = --, (51) 

7 7 



in which we recover 



^'-^^^(IT?' (52) 



as expected. 



S. 7. Multiple absorbers 

If we have more than one absorbing element at positions X\, x% . . ., we will use the same 
formula for the absorption efficiency (150"]) . but with transmission matrix 

T toiai = m{x 2 - x l )T 2 $(x 3 -x 2 )--- T N . (53) 

Here, Tj describes the absorption properties of a given absorber and thus depend on the 
value of 7i, while the phase matrix $(x) is given by 



e 



£00) Q 



^ = e -*« (54) 
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Figure 3. Detection efficiency. (Top) Absorption probability vs. effective decay rate 
for a setup with N = 1,2,4 and 8 qubits (black, green, blue, red) either in cluster 
(dashed) or array (solid) regime. In the first case, where qubits are close together, the 
efficiency is fundamentally limited to 50%, while in the other case there is no upper 
limit. In both cases, having more qubits allow us to use weaker couplings or stronger 
decay rates. (Bottom) Detection efficiency versus the separation d of the absorbing 
elements for an array of periodically distributed Josephson junctions. All curves have 
been produced for resonant qubits. 

and contains the accumulated phase of the electromagnetic wave when travelling between 
consecutive absorbers 

8(x) = —x = —x. (55) 
A v 

2.8. Analysis of the efficiency 

The relation ( 1471) is valid for one, two and more absorbing elements, placed at different 
positions. If we have a single element, we have already found the absorbed fraction (|52|) . 
As shown in Fig. [3l this quantity is upper bounded by 50% which is reached for resonant 
qubits with 7 = 1. This means that a single absorbing element, without knowing the 
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Figure 4. Dependency of the detector efficiency on the parameter 7 and on the 
separation between absorbers, in units of the microwave wavelength, d/X. We have 
chosen N — 4 (above) and N = 16 (below). Note that the maximum gets broader with 
N, but it is always achieved at d — A/4. The solid lines denote the cuts corresponding 
to the plots in Fig. [3] 



arrival time and no particular external modulation, can only detect in average half of 
the incoming photons. 

The situation does not improve if we place N absorbing elements close together 
[Fig. [lb] . When the wavelength of the microwave field is large compared to the size of 
the absorber cluster, the total transfer matrix is given by 

T = T»=(i + -a) =1 + -A. (56) 

V 7 / 7 

This is just the same matrix with a smaller effective decay rate 7/iV, and the total 
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efficiency remains limited to 50%, 

aN= [l + 7 /iVp -2- (57) 
We can do much better, though, if we place a few absorbing elements separated by 
some distance d [Fig. [lb]. In this case the total transfer matrix is given by 



T N {d) 








(5f 



The total efficiency now depends on two variables, 7 and the phase 9 = Uod/v or the 
separation between absorbers, d. As Fig. 0] shows, the optimal value of the phase is 
6 = 7r/2, corresponding to A/4 separation and, in this case, the total efficiency is no 
longer limited. For instance, as shown in Fig. [3] for two and three qubits on-resonance 
the efficiency can reach 80% and 90%, respectively. Not only it grows, but it does so 
pretty fast. 

An interesting feature is that, as we increase the number of qubits, the absorbed 
fraction becomes less sensitive to the qubit separation, which allows for more compact 
setups than one would otherwise expect. For instance, in Fig. [5] we show the total 
size of a setup computed for a fixed detection efficiency and a given number of qubits. 
Since the optimal separation behaves as d ~ 1/N, the total system size for a fixed 
detector efficiency, 78%, 80% and so on, remains bounded, even for large number of 
qubits. Moreover, as one increases the number of qubits, the efficiency grows rapidly. 
We can define the value 

aoo(£) = lim maxa(£, N, 7, d = C/N), (59) 

N^oo 7 

which gives an idea of what is the maximum detection efficiency for a given circuit size. 
The value shown in Fig. [5b approaches the limit of 100% quite fast and gives us an idea 
of the minimal size of a detector which is needed to obtain a given efficiency. 



2.9. Robustness against imperfections 

There are many factors that will condition the actual efficiency of a photodetector. Some 
of them will have to be discussed later on in the context of the proposed implementation, 
but others can be analyzed already with the present theory. 

The first source of errors that one may consider are systematic differences in the 
fabrication and tuning of the absorbing elements. These fluctuations are currently 
unavoidable, and may even evolve through the lifetime of a setup, due to changes in the 
temperature, fluctuations of impurities, among others. These systematic errors could be 
modeled by random perturbations in the parameters of the three-level systems, either 
due to inhomogeneous broadening (different frequencies u>i), inhomogeneous decay rates 
(different Tj) or changes in the coupling strengths (Vi). However, since the scattering of 
photons is described by a single parameter 7^ per qubit, it is more convenient to model 
the errors as random changes in these values. 




Figure 5. Optimal separation of detectors, (a) Asymptotic detector efficiency a^L) 
as a function of the setup size L, and computed as in Eq. ([59]). (b) Minimal size L of 
a setup with N qubits required to attain a detection efficiency a = 78% (dash-dot), 
85% (dashed) or a = 92% (solid). 



We have performed a numerical study using a constant probability distribution, 
P(^y) = l/2<7 7 7 for I7 — 7I < <t 7 7, with a maximum relative error <j 1 that ranges from 
up to 20%. Using 10000 random realizations with different parameters {jijiLi, we 
made statistics of the detector efficiency, a. The results are shown in Fig. [61 where we 
plot the standard deviation of the maximum detection efficiency, a Q , given by 

*« = V / ((« 2 )) " ««» 2 , (60) 
where the double bracket ((•••)) denotes average over the random instances. According 
to our results, the maximum efficiencies for N = 4,8, 16, and 32 are a = 
92.4%, 97%, 98.9% and 99%, and in the worst case of 20% maximum systematic errors, 
where one standard deviation implies an efficiency decrease of 0.7% for 4 absorbers and 
0.1% for 32. In other words, we find up to two orders of magnitude smaller errors in 
the detector than in the fabrication. 

2.10. Detector bandwidth and dephasing 

There are two other sources of error which are intimately connected. One is the 
bandwidth of the incoming photons, Au, which can be related to the length of the 
photon wavepacket, r ~ 1/Au. Since the detector efficiency depends on 7 oc T+i(u— uq), 
which contains the detuning between the photon and qubit frequencies, some spectral 
components may only be weakly detected. To solve this problem we need the detector 
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Figure 6. Sensitivity of the detector to systematic errors in the absorbers. We plot the 
standard deviation of the absorption efficiency, a a , in percentiles, vs. the maximum 
systematic random error, <r 7 , averaged over 10000 random samples with N = 4, 8, 16 
and 32 qubits (top to bottom). 

efficiency to be weakly sensitive to variations in the detuning of the photon, 5 = uj — u> , 
that is we need a broadband detector. 

To estimate the bandwidth of the photodetector we have applied the following 
criterion. We compute the optimal value of 7 op4 for which the detection happens at 
maximum efficiency with zero detuning, 5 = uo — uj = 0. We then look for the two 
values of the detuning at which 7 = 7 op t(l + U5/T) causes a reduction of the efficiency 
of a given percentage. As Fig. [7b shows, the bandwidth grows with a power law N u , 
with an exponent which is unfortunately not too large. However, Fig. [7^ shows that 
the detuning may actually increase the efficiency, probably indicating that the previous 
analysis is too limited and that the detector design may involve optimizing both T and 
S not only for achieving a certain efficiency, but also to increase the efficiency. 

Another related problem is dephasing. As we will discuss later on, quantum 
circuits are affected by 1/f noise. Part of this noise can be understood as oscillating 
impurities that change the electromagnetic environment of the absorbers, and thus 
the relative energies of the and 1 states. This error source is modeled with a term 
ea z = e(|l) (1| — |0) (0|), where e is a random variable — either classical, or quantum, 
from a coupling with the environment — . When one averages over the different noise 
realizations, the result is decoherence. 

We can get rid of the noise for each realization using a unitary operator U(t) ~ 
exp(— i Jq ea z dr/h). When we use this operator to simplify the Schrodinger equation, 
creating the equivalent of an interaction picture, the result is that we can translate the 
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Figure 7. (a) Change of the detector efficiency as we introduce a detuning between the 
photon and absorber frequencies, (b) Bandwidth of the photo detector vs. the number 
of absorbers. Given the optimal working point j pt which maximizes the detection 
efficiency at zero detuning, 6 = lo — lu , we look for the value of the detuning such that 
the efficiency at 7 = 7 op t(l + is reduced by 5% (circles) or 20% (crosses). The 
solid lines are fits to iV 016 and N - 2 . 



accumulated random phase, 



edr/h, 



0(t,e) = 
into the coupling 

H mt ~ V I 5(x) 



A + e^ipi |l)(0| + H.c. 



(61) 



(62) 



In other words, from the point of view of the absorber it is as if the incoming photon was 
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affected by a random, but slowly changing phase [j] that may either shift the frequency of 
the photons or broaden their spectral distribution. In either case, if the photodetector 
has a large enough bandwidth we may expect just a minor change in the detector 
efficiency 

3. Model implementation: microwave guide with phase qubits 

In this section we detail a possible implementation of our scheme which is based on 
elementary circuits found in today's experiments with superconducting qubits. We will 
show how our previous theory relates to the mesoscopic physics of these circuits and 
compute expressions for the relevants parameters, v, V, j, . . . in terms of the properties 
of these circuits. 

For the waveguide we will consider a coaxial planar microwave guide such as the 
ones employed to manipulate and couple different qubits [T9J HI] and described in detail 
in Sec. HJ Note, however, that unlike in Ref. [19] our waveguide will not be cut at 
the borders and it will not form a resonating cavity. For the bistable elements we 
will consider a superconducting qubit, the so called "phase qubit" or "current-biased 
Josephson junction" (CBJJ), which has a set of metastable levels that, by absorption 
and emission of photons, may decay to a different, macroscopically detectable current 
state. 

3.1. Current-biased Josephson junction 

As mentioned before, our detection element will be a CBJJ. The model for this circuit is 
shown in Fig. [HJ there is a Josephson junction shunted by an current source which can 
be modeled by a very big impedance. The bias current / causes a tilting of the energy 
potential in the junction, creating metastable regions with a finite number of energy 
levels, that tunnel quantum-mechanically outside the barrier. 

The quantization of this circuit renders a simple Hamiltonian [12] 

H=^-Q 2 + U( ( j ) ), (63) 

expressed in terms of the charge in the junction and the flux 0, at a node of the circuit 
[Fig. [HJ. The Hamiltonian contains the usual capacitive energy, expressed in terms of 
the large capacitance of the junction, Cj, and a potential energy due to the inductive 
elements. Modeling the current source as a large inductor with a total flux that supplies 
a constant current, = J, we obtain a highly anharmonic potential 

U{4>) = -I <p cos(0/^ o ) - H- (64) 

Note that even though the actual flux quantum is 0o — h/2e, in order to avoid 2n factors 
everywhere it is convenient to work with ip = o /2vr. 

The quantization of this model corresponds to imposing the usual commutation 
relations between the canonically conjugate variables, the flux <p and the charge Q, that 

| Remember that the noise source is 1 /f and dominated by low frequencies 
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Figure 8. (left) Scheme for a current-biased Josephson junction, (right) Energy levels. 
There are N s metastable energy levels with anharmonic frequencies. 



is [0, Q\ = ih. Given the relevance of the anharmonic terms in Eq. (|64l) it soon becomes 
evident the convenience of working in the number-phase representation <p = <£>n x 8, and 
Q = 2e x N with operators that satisfy [8, N] = i. In these variables the Hamiltonian 
becomes 



H = E C N 2 - I <po 

with the junction charging energy 
(2e) 2 h 2 



cos(fl) + —8 



E, 



c 



2C, 



(65) 



(66) 



2<f%Cj 

3.2. Harmonic approximation 

When the bias current I is very close to the critical current Jo, we have the situation in 
Fig. [HI in which the junction develops a metastable, local minimum of the potential at 
8 close to 7r/2. It is then customary to approximate the potential by a cubic polynomial 
and describe the dynamics semiclassically, with a coherent component that describes the 
short-time oscillations around the local minimum and a decay rate to the continuum of 
charge states which are outside this unstable minimum. 

The semiclassical limit is characterized by just two numbers, the plasma frequency 
of the phase oscillations around the minimum 

1/4 



1- 4- 



and the barrier height 



AU 




(67) 



(68) 



that prevents tunnneling outside this minimum. Using semiclassical methods it is 
possible to estimate the number ~ N s = AU/huj p of metastable states in this local 
minimum, and approximate their energy levels, 

.anh 



E n /h = nojp + oo° 



iT r 



(69) 
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including an anharmonic component, u^ nh , and an imaginary part which is the rate at 
which the state decays into the continuum 

Tn=ui m^ e ^ {70) 



In this work, we use the fact that for two consecutive levels these rates can be very 
different, r n+1 /T n ~ 1000. If the decay rate of the first level I\ is large enough while 
still keeping T small, we can treat the levels |0) and |1) as the two levels of our absorbing 
element, the continuum playing the role of the detectable states \g) . With this we have 
now the parameters uj = uj p and r = F 1 of our model. It only remains to find out the 
coupling V and the group velocity v. 

In order to quantify the coupling between the junction and the microwave guide we 
will need to explicitely make the calculations leading to the previous results and find 
out the harmonic approximations for the lowest energy levels. We begin by noticing 
that the energy minima are reached at a value of the phase 

sm(6 m ) = -f =: r (71) 

which is very close to n/2, as predicted. Around this minimum, 8 = 9 m + 6, a harmonic 
approximation gives 

H ~ E C N 2 + \u"6 2 , (72) 
with a curvature of the potential 



U"(6 m ) = JoV?o cos(fl m ) = Jo^oVl^. (73) 
To diagonalize this Hamiltonian we introduce operators 

# = -Li(ot - o), (74) 

=^(a + at), (75) 
with usual commutation relations, [a,a^] = 1, and impose 

huo p . (76) 



2a 2 4 2 
This provides us with the plasma frequency, 



but also with the "strength" of the number and phase operators, 



2 2E c E c . . 

a= iu^ = ^ (78) 

which is related to the charging energy and the plasma frequency in a simple way. 
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3.3. Qubit-line in the dipole approximation 

Given that the total capacitance of the transmission line is much larger than that of 
a single superconducting qubit, we can apply the dipole limit to study the coupling 
between both elements and assume that the interaction term goes as 

H = q-^—V(x,t). (79) 

Following the theory in Sec. HJ the potential created inside the waveguide can be written 
as a function of the charge distribution on the conductor 

(80) 

c ox 

Using the eigenvalue equation d x Wk{x,t) = ikwk{x,t) and the dispertion relation (|89|) . 
we obtain 




V(x, t) =iJ2\l -^7 a k w k (x, t) - alw k (x, t)* 



11) 



If we assume that the incident and outgoing wavepackets have a small bandwidth, 
uj k ~ ujq, we obtain the potential as a function of the propagating fields 



V 



, v ^r+^-lPl-tPj]. (82) 

As a consistency check, we can consider the case of a small transmission line, forming 
a "cavity" or resonantor in which the modes are very well separated. In that limit we 
only need to consider a single momentum, p, and everything reduces to the formula 

V(x,0)=t^(a p -al). (83) 

Here, C = L x c is the total capacity of the transmission line and the quantity 
(hu!k/2C) 1 / 2 is the r.m.s. voltage. Note that, consistently with the expression in the 
experiment of Wallraff et al. [HI [19] , the coupling between the qubit and a particular 
mode is inversely proportional to the square root of the microwave guide length, L. 
However, the coupling constant with the fields in Eq. (|82|) is insensitive to the total size 
of the transmission line. 



3.4- Final parameterization of the setup 

Following the previous considerations we will write the interaction Hamiltonian between 
the junction and the microwave field as in Eq. (1791 using the charge operator fTM|) . The 
coupling between the stripline and the two lowest levels of the CBJJ is proportional to 
the constant 



V=^ r iJ^. (84) 

Let us remind that u is the photon frequency, C g is the gate capacitance between the 
junction and the waveguide, c is the capacitance per unit length of the waveguide and 
a is the harmonic oscillator wavepacket size (1781) . 
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In addition to the coupling strength we need to compute the group velocity v. 
The precise derivation is left for Sec. HI but we advance that this can be written using 
properties of the waveguide 
1 1 



. (85) 
Vd cZ 

where I is the inductance per unit length and Z the impedance. When introduced in 
our scattering problem, this results in the effective decay rate 7 

- —hr - n / cZ o hr _ a 2 h T 
7 ~ V 2 c\ 2 e 2 hu /a 2 c c? 2 e 2 Z, cu ' 1 j 

expressed in terms of ci 2 = C g /(C g + Cj) for resonant qubits (i.e. no detuning). We 
want to remark that there are plenty of parameters to tune, so that realistic conditions 
of high efficiency are achievable. 

In practice one would begin by fixing the number of superconducting absorbers 
that can be used, N. This fixes a value of 7 for which absorption is optimal and which 
determines the optimal decay rate of the first energy level 

r = f^o- (87) 

As an example, let us compute the optimal configuration for one qubit. Using the 
numbers in Berkley's experiment [15] . we will fix the qubit capacitance Cj = 4.8pF 
and C12 = 0.13. The charging energy is about Ec/h = 0.10 GHz, which together with a 
plasma frequency uj = 5 GHz gives a 2 = 0.02. Putting the numbers together 

r = 5i!?r°- < 88) 

so that for inductances from 10Q to 100Q this gives T = 10 — 100 MHz. Higher decay 
rates can be achieved by increasing C g or decreasing the charging energy. For instance, a 
factor 2 increase of C g gives a factor 4 increase in T which now ranges from 30-300 MHz. 



4. Photons in transmission lines 



In this section we study the mathematical models and physical properties of a coaxial 
coplanar microwave transmission line. This device consists on a central conductor 
carrying the wave and enclosed by two conductors set to the ground plane. A simple 
but effective model for such a line is obtained by coupling inductances and capacitors as 
shown in Fig. [9j As we will show later on, if we denote by I and c the inductances and 
capacitances per unit length, this model predicts the propagation of microwave fields in 
the transmission line according to a dispersion relation 

Uk = v\k\ = —p=\k\, (89) 
vet 

where fik is the momentum of the photon, Uk is the frequency and v is the group velocity. 
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c, 



C r 



Figure 9. A transmission line can be modeled as a set of inductances and capacitances. 
4-1. Discrete model 

In order to analyze this circuit from Fig. [9j we must write down the Kirchhoff 's law for 
the n-th block containing 4 nodes. When combining all equations and leaving as only 
variables the branch intensities, we obtain the set of second order differential equations 



-L - 



d 2 L 



-[2/ n 



'n-l 



'n+lj 



(90) 



dt 2 C ' 

These equations are similar to those describing an infinite set of oscillators of mass 
m oc L and spring constant k oc 1/Cq. In analogy with the mechanical case, if we 
assume periodic boundary conditions to better reproduce propagation of charge, we 
find travelling wave solutions 

2tt 



I n (t)=I e** x »- ut >, k= — xZ, (91) 

where x n = a x n, a is a parameter denoting the distance between neighbor oscillators, 
p is the momentum of the wave and L is the length of the line. A direct substitution of 
this expression in Eq. ()9Up . gives the dispertion relation 

i 1 

2 



u(k) 



[1 — cos(fca)) 



LqCq 



k 



(92) 



which is approximately linear for small momenta, long waveguides or thin discretization. 
Using the inductance and capacitance per unit length 

I = L /a, c = Co/a, (93) 

we obtain the group velocity v and dispersion relations introduced before (1591 . 



4-2. Lagrangian formalism and continuum limit 

The previous evolution equations (19"01 can be obtained from the Lagrangian 



L 



E(y« 



1 



2C 



(Qn ~ Qn+l) 



using the Euler-Lagrange equations 
d ( dL\ dL 



dt \dQ n ) dQ n 



(94) 



(95) 
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A more realistic model for a continous transmission line is obtained by taking the 
limit of infinitesimally small capacitors and inductors, a — > 0, while preserving the 
intensive values I and c. In the continuum limit we replace the discrete charges by 
a continuous charge density distribution, Q n ~ a x q(x n ,t), which should be a smooth 
and integrable function of the position variable x. Through this procedure one obtains 
a total Lagrangian expressed in integral form 



dx 



l -{d t qf - ±-(d x qf 



(96) 



The charge distribution described through the field q(x, t) obeys a wave equation 

%q(x,t)-\%q(x,t) = 0, (97) 
where the group velocity is precisely the one already found in the discrete case. 

4-3. Hamiltonian and quantization 

The quantization of an electrical circuit is a three-step process. We begin by identifying 
the classical canonically conjugate variables: the charge distribution q(x,t) and the 
associated momentum 
dC 

U q = — = lq(x,t). (98) 
The Hamiltonian is then built using the prescription 

\2 (Si „\2~ 

(99) 



H — J dxll q q — L = J dx 



21 2c 



Finally, the variables {q,U q }, are replaced with Hermitian operators {q, tl q }, with 
commutation relations 

[q{x, t), fl q (x', t)} = ih8{x - x'). (100) 

The previous Hamiltonian is quadratic and it can be diagonalized. We begin by 
analyzing the dynamics of these operators in the Heisenberg picture and noticing that 
the charge satisfies a wave equation (|97|) . It therefore makes sense to look for a set of 
normal modes made of plane waves 

q{x,t) = J2N k q k (t)u k (x), (101) 

k 

1 2n 
u k {x,t) = —=exp(ikx), k G — x Z, (102) 
v C £ 

where q k are our new dynamical variables and N k is a normalization constant which will 
be fixed later on. Note that we are using periodic boundary conditions because they are 
best suited for describing transport, but we have not fixed the transmission line length 
C, which may be arbitrarily large. 

Replacing the expansion above in (J97j) . produces second order differential equation 
for the unknown operators 

'q k (t)+oo 2 k q k (t) = 0. (103) 
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Since both q and q k are physical observables and Hermitian operators, the most general 
solution is q k (t) = a k e~ luJkt + a\,e luJkt , expressed in terms of a time independent operator, 
a k . The commutation relation f llOOp produces both bosonic commutation relations for 
the a k 

[a k ,a\i\ = S k k>, [ak,a k '} = [%,%,] = 0, (104) 

and N k = ft/(2u k l). Finally, with the orthonormalized wave functions w k (x,t) = 
Uk(x) exp(—iu k t), we arrive at 



q(x, t) =Y / J^- l [a k w k ( Xl t) + H.c] (105) 



ftqOM) = i \l~Y~~ [ a l w k(x,ty - H.C. 

Given the specific form of the canonical operators, we may obtain a particular 
dispersion relation uo k that diagonalizes the Hamiltonian. Using the relations 

w k {x,t)w k ,(x,t)dx = We- i( " fc+ " fc,)t , (106) 

w k (x,t)w k >(x,t)*dx = S k -k>, (107) 
and imposing 

Hud 1 ft 1 , 9 / \ 

— - x — = x — x k , 108 

2 21 2u k l 2c y ' 

we will be able to cancel all terms proportional to a k a- k and a\a}_ k , obtaining a set of 
uncoupled oscillators 

H = J2ftUk(ala k + ~y (109) 
where the dispersion relation is strictly the one introduced before in Eq. (1591 . 

4-4- Linearization 

In our work we focus on states that contain photons with momenta around \ko\ or Wq/v, 
where ujq is the principal frequency of the wavepacket. We thus introduce two field 
operators representing the right- and leftward propagating photons, 

ip r (x,t) = Y a k w k (x,t), (110) 
keB 

?pi(x,t) = ^2a_ k w_ k (x,t), (111) 

keB 

where B = [k — A, k + A] is the desired neighborhood around the principal momentum, 
characterized by a sensible cut-off A. These two fields satisfy the evolution equations 

id t ip r (x,t) = — iftvd x ip r (x,t), (112) 
id t ipi(x,t) = +iftvd x i[)i(x,t), (113) 
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and we can write the effective Hamiltonian 

H= ipl(—ihvd x )ip r + ipii+ihvd^ipi dx. 



(114) 



There is a catch in the previous expansion: the operators ijj r j do not have the usual 
bosonic commutation relations. For instance 



where the right hand side represents a truncation of the distribution 8{x — y) within 
the interval of momental B. In this work we are interested in the asymptotics of the 
absorption process, treating long times and long enough wavepackets for which the 
rotating wave approximation is valid. In this case we are justified to approximate 



and treat the right and left propagating fields as causal. 
5. Conclusions 

We have developed the theory of a possible microwave photon detector in circuit 
QED, and studied diverse regimes, advantages and difficulties with a realistic scope. 
Though we believe that our contribution will boost the theoretical interest and 
possible implementations in microwave photodetecion, we will summarize the potential 
limitations and imperfections of our proposal. First, the bandwidth of the detected 
photons has to be small compared to the time required to absorb a photon, 
roughly proportional to 1/T. Second, the efficiency might be limited by errors in the 
discrimination of the state \g) but these effects are currently negligible [29]. Third, 
dark counts due to the decay of the state |0) can be corrected by calibrating T 
and postprocessing the measurement statistics. Fourth, fluctuations in the relative 
energies of states |0) and |1) , also called dephasing, are mathematically equivalent to 
an enlargement of the incoming signal bandwidth by a few megahertz and should be 
taken into account in the choice of parameters. Finally, and most important, unknown 
many-body effects cause the non-radiative decay process 1 — >• 0, which may manifest 
in the loss of photons while they are being absorbed. In current experiments [29], this 
happens with a rate of a few megahertzs, so that it would only affect long wavepackets. 

Our design can be naturally extended to implement a photon counter using a 
number of detectors large enough to capture all incoming photons. Furthermore, our 
proposal can be generalized to other level schemes and quantum circuits that can absorb 
photons and irreversible decay into long lived and easily detectable states. 

We expect to have contributed to the emerging field of detection of travelling 
photons. Its success may open the doors to the arrival of "all-optical" quantum 
information processing with propagating quantum microwaves. 




(115) 




(116) 
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